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Abstract 

^ ■ The problem of scattering of linear internal waves from small com- 

5/3 \ pact sea bottom inhomogeneities is considered from the point of view 

of mode-to-mode scattering. A simple formula for modal conversion 
coefficients C nm , quantifying the amount of energy that is scattered 
r- 1 , into the n-th mode from the incident field m-th mode, is derived. In 

Q_ij this formula the representation of inhomogeneities by their expansions 

into the Fourier and Fourier-Bessel series with respect to angular and 
radial coordinates respectively are used. Results of calculations, per- 
formed in a simple model case, are presented. The obtained formula 
■ can be used for a formulation of the inverse problem, as it was done 

CN . in the acoustic case [21 [3] . 

| Keywords: internal wave, scattering 

o 

1 Introduction 

o 

The concept of mode-to-mode scattering was considered in the context of 
! the acoustic scattering from small compact irregularities of the ocean floor by 

Wetton, Fawcett [T]. In their work some simple formulas for modal conversion 
coefficients, quantifying the amount of energy that is scattered from one 
^ normal mode of the sound field to another, were derived. Recently new 

formulas for these coefficients were obtained by Zakharenko |2] and applied 
to the inverse scattering problem in the subsequent work This paper 
contains the detailed derivation of such formulas in the case of scattering of 
linear internal waves from small compact sea bottom inhomogeneities. Some 
numerical examples are presented. 
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2 Formulation and derivation of the main re- 
sult 

We shall use the linearized equations for inviscid, incompressible stably strat- 
ified fluid, written for the harmonic dependence on the time with the factor 

e ~iuit j n f orm 

-iup u + (3P X = 0, 
-iup v + (3P y = 0, 
-iup w + (3P z + (3 P i = 0, (1) 
-iupi + wp 0z = 0, 

U X + Vy + W z = 0, 

where x, y, and z are the Cartesian co-ordinates with the z-axis directed 
upward, po = Po(z) is he undisturbed density, p\ = pi(x,y,z) is the pertur- 
bation of density due to motion, P is the pressure, and u, v w are the x, 
y and z components of velocity respectively. The variables are nondimen- 
sional, based on a length scale h(a typical vertical dimension), a time scale 
N^ 1 (where N is a typical value of the Brunt- Vaisala frequency), and a den- 
sity scale p (a typical value of the density). The parameter f3 is g/(hN 2 ), 
where g is the gravity acceleration. 

The boundary conditions for these equations are 

w = at z = fl, 

(2) 

w = —uH x — vH y at z = —H, 

where H = H(x, y) is the bottom topography. 

We introduce a small parameter e, and postulate that the components of 
velocity and the pressure are represented in the form 

u = u + eui + . . . , v — v + evi + . . . , 

w = w + ew\ + . . . , P = Pq + ePi + . . . 

We suppose also that the bottom topography is represented in the form 
H = ho + ehi, where ho is constant and h\ = hi(x,y) is a function of x, 
y vanishing outside the bounded domain Q, which in the sequel is called a 
domain of inhomogeneity. 

Excluding from the system ([T]) p 1 and substituting the introduced expan- 
sions, we obtain 
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-iup (u + eiti + . . .) + (3{P 0x + eP lx + ...) = 0, 
-iup (v + evi + ...) + (3(P 0y + eP ly + ...) = 0, 
(u 2 p + (3p 0z )(w + ew x + ...) + iu(3(P 0z + eP zl + ...) = 0, ° ] 
(u 0x + eu Xx + •••) + ( v o y + e^iy + ...)+ w 0z + ew Xz + . . . = 0, 
with the boundary conditions 

w + ew\ + . . . = at z — 0, 
w + ewi + . . . = -Oo + eu x + . . .)(h x0 + eh lx ) (4) 

—e(v + evi + ...)(hiy + ...) at z = -H. 

Separating terms in various orders of e, we obtain a sequence of boundary 
problems. 

At order O(e ) we have 

-iujp u + (3P 0x = 0, 

-iujp v + (3P 0y = 0, ^ 
(u; 2 po + f3po z )wo + iuPPoz = 0, 
UOx + Vo y + w 0z = 0, 

with the boundary conditions 

w = at z = 
wo = at z = —h . 

Differentiating the third equation in §5§ twice with respect to x and twice 
with respect to y, summing obtained equations and replacing /3(Po zxx + Po zyy ) 
by — iu(poWo z ) z , we obtain 

(uj 2 p + pp 0z )(w 0xx + w 0yy ) + uj 2 (p w 0z ) z = 0. (6) 

We seek a solution to this equation in the form of the sum of normal modes 
Wq = e l ( kx+lyS) (j)(z), where <p is the eigenfunction of the spectral boundary 
problem 

-{u 2 p + (3p 0z ){k 2 + l 2 )<p + u 2 { Po <f) z ) z = 0, 

0(0) = <f>(-ho) = 0, U 

with the eigenvalue A = k 2 + I 2 . It is well known that the problem ([7]) has 
countably many eigenvalues A„, which are all positive. The corresponding 
real eigenfunctions <p n we normalize by the condition 

(uj 2 p + pp 0z ) (f) 2 dz = / p o (0 2 ) 2 dz = 1. (8) 

-h K + I J_ hQ 
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The eigenfunctions 4> n and <p m with n m are also orthogonal 

(0n,0m)=O (9) 

with respect to the inner product 

r° 

(0, 1>) = - (^Po + Ppoz) H> dz (10) 

J— ho 

In our scattering problem wq is the incident field, and we shall calculate the 
main term of scattering field W\, so we act in the framework of the Born 
approximation. 

At the first order of e we obtain the following system of equations: 

-iupoUt + pP lx = 0, 
-iup v 1 + (3P ly = 0, 
(u 2 p + (3p 0z )w 1 + iu(3P lz = 0, 
Ui x + v ly + w lz = 0, 

with the boundary conditions 



wi = at z = 0, 
Wi = —Uohi x — Vohiy at z = —h — eh\. 



(12) 



So far as we are interesting in the connection of modal contents of incident 
and scattering fields, we suppose that the incident field consists of one mode 
Wq = e l ( knX+lny ' > (j) n (z) . Reducing the second boundary condition (12) to the 
boundary z = —h with taking into account the explicit form of Wq, we obtain 
the new boundary condition for Wi at the boundary z = —h : 



Wi 



lk n tin , \ i(k n x+l n y) i , ;, 



Reducing the system fill I) in the same manner as it was done for the 
system (jSJ), we obtain the equation for W\. 

(u 2 p + (3p z)(wixx + Wiyy) + uj 2 (p w lz ) z = 0. (14) 

We seek the scattering field in the form w\ = J2m=i C n m{x-,y)4>m-, the func- 
tions C nm (x,y) are called the modal conversion coefficients. To obtain the 
equation for C nm we substitute the postulated form of W\ to the ffT3"]) . multipli- 
cate it by the function <fi m and integrate from — h to 0. Using the conditions 
of orthogonality and normalization (JHJ), (J9J) and the boundary condition ( TTSl . 
we finally obtain 
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q2 q2 

where 

Writing the solution to the equation (TT5]) as the convolution of the 
fundamental solution (Green function) of the Helmholtz operator G = 
(—i/4)HQ 1 \y/k% l + l^R) with the right-hand side F, we have 

C nm (x r ,y r ) = -\f f FH^(y/kl + P m R)dydx, (16) 

x y 

where R = \J(x — x r ) 2 + (y — y r ) 2 and by the index r we designate the point 
of registration of the field. 

Integrating by parts the terms containing h\ x ,hi y and passing to the 
cylindrical coordinate system with the origin in our domain of inhomogeneity 
and such that k n = K n , l n — 0, x — rcosa, y = rsina, we obtain 

oo 2tt 

l«m ~ f f . v .,(1) 



C nm = --^G / /i ie lK " rcosQ cos(^-a r )^ iJ (fi: mJ R)rrfarfr, (17) 
4k„, 







G = uJ 2 po4>nz(-h )(f) mz (-h ), R = A/r 2 + r 2 - 2rr r cos(a - a r ), (ry, a r ) are 
the polar coordinates of the registration point, tan(^) = r sin(a — a r )/(r r — 
r cos(a — a r )). 

Using the addition theorem for the Bessel functions we express contained 
in (TT7]) cos ifiHi (K m R) and sin ipHi (K m R) in the form: 

oo 

{sinC0) } ^1 \ K mR) = ^2 H kll( K rnrr)Jk{K m r) ^y s i n k( a -aj) } ' 

fc=— oo 

From now on we shall assume that the distance r r to the registration point 
is big enough to replace the functions if£_i(/c m r r) by their asymptotics 

H^l x {K m r r ) ps v / 2/(vr/€ m r r )exp [i(/« m r r - (ir/2)(k + 1) - tt/4)] . 

Then, expanding /ii(r, a) as function of a in Fourier series with the coeffi- 
cients hi u (r), after integration with respect to a, we obtain 
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C nm = !^ ^«p(^v - W4) GcosQr £ ^ e .^ 



oo „ 

)Ju+k(K n r)r dr 

k=— oo q 

Changing the order of integration and summation we can achieve further 
simplification by using the formula 



Y J k {n m r)Ju+k{K n r)e lkar = J v {£,r)e 



ivB 
k=— oo 



where £ = \/ x? m + — 2x m x„ cos a r , 6 = arctan — — — — . We ex- 

t\"Yi f^m COS OLf 

pand now the radial coefficients h\ v (r) on the segment [0, L], where they do 
not vanish, in the Fourier-Bessel series 



Mr) = J, (f r) 
p=i ^ ' 



where 7^ are the positive roots of the function J u , J v (jp) = 0. Substituting 
this expansion in (1181) and taking into account that 



r 2 - ' 

we obtain the final expression for modal conversion coefficients 







iL 2 V2n y/K^exp(iK m r r - in/ 4) 

^nrn = ~j= Cr COS Ct r 



(19) 



WUHL)e-^^Y.fv 



X 

v=— CO 



P =l [ P ^ 
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3 Numerical examples 



For a model example we choose p = e Xz ,f3 = A 1 and H = 1 . Then the 
spectral boundary problem is written in the form 

U 2 <fizz — U 2 \(f)z — ^(id 2 — 1)0 = 0, 

0(0) =0, 0(-l) = O. 

The eigenfunctions of such a problem are = Ae Xz l 2 sin((/ + l)nz) with the 
eigenvalues 

+ 1)% 2 + A 2 /4 



K 



Here A = y/2/ — a; 2 ) by the condition (9). For the calculations the value 
of parameter A was taken to be equal to 0.003 , which corresponds to the 
typical stratification in the ocean shelf zones. The domain of inhomogeneity 
has the form of the ellipse with the big and small radii a and b of which were 
taken in proportion a : b = 2 : 1, and in this region 



hfay) = 0.05^1- 

In the figure are presented the results of calculations with uo = 0.5 and the 
angle of incident field ao = 0, conducted for various wave sizes kcl of the 
scatterer. We note that according to the meaning of the small parameter 
e, in these calculations e = 0.05. For the presentation of results we use the 
scattering amplitude 
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Figure 1: Absolute value of scattering amplitude: na — 1 (a,b), na = 2 (c,d), 
rea = 8 (f,g) 
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